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ABSTRACT 

Infrared spectra from the Spitzer Space Telescope (SSC) of many debris disks are 
well fit with a single black body temperature which suggest clearings within the disk. 
We assume that inside the clearing orbital instability due to planets removes dust 
generating planetesimal belts and dust generated by the outer disk that is scattered 
or drifts into the clearing. From numerical integrations we estimate a minimum planet 
spacing required for orbital instability (and so planetesimal and dust removal) as 
a function of system age and planet mass. We estimate that a 10 8 year old debris 
disk with a dust disk edge at a radius of 50 AU hosted by an A star must contain 
approximately 5 Neptune mass planets between the clearing radius and the iceline 
in order to remove all primordial objects within it. We infer that known debris disk 
systems contain at least a fifth of a Jupiter mass in massive planets. The number of 
planets and spacing required is insensitive to the assumed planet mass. However an 
order of magnitude higher total mass in planets could reside in these systems if the 
planets are more massive. 



1 INTRODUCTION 



Recent observations made with the Spitzer Space Telescope 



circumstellar disks dRieke et al. 2005; Bcichman et al.ll2005l; 


Gorlova et al.1 2006: 


Chen et al.l [20061; iBrvden et al.1 [20061; 


Beichman et al. 2006 


). Theses disks have total opacity, as 



estimated from the fraction of stellar light re-emitted in the 
infrared, in the range 10 -3 - 10 -5 . The fraction of stars with 
disks increases from 2% for r > 10" 4 to 12% for r > 10~ 5 
([Brvden et al.ll2006l ). and the fraction of stars with disks 
detec t ed increases with de creasing stellar age l|Rieke et al.l 
120051 ; [ Gorlova et al.l 120061 ) . Infrared spectra of many de- 
bris disks are well fit with a single black bo dy tempera- 



ture suggesting that they possess inner holes ( Chen et al.l 
2006). Planets are expected to be responsible for the inner 
clearings, however the number and masses of these planets 
has yet to be constrained by dynamical arguments. An ex- 
amples of a nearby system with an inner clearing that has 
been resolved wit h imaging is the Fomalhaut system (e.g., 
iKalas et al.ll2005l ). 

With the discovery of multiple planet extrasolar 
systems, the stability of planetary systems has been 
inves tig ated with renewed intere st (e.g., 



Ford et al 



200l| ; iBarnes fc Raymond! 12004 iLepage & Duncan 



2004; Ford et al 



2005 



Barnes fc Greenberd 120061 : 



Raymond fc Barnes! 2005[ 
IChatteriee et al.1 120070 



A 

system of two planets on low-inclination and eccentricity 
orbits will never experience mutual close encounters if the 
initial semi-major axis d ifference, A m , m easured in mutual 
Hill radii , exceeds 2\/3 jGladmanlll993T ). For multi-planet 
~aU 



systems, IChambers et 



(|l996l ) 



numerically investigated 



the dependence of the timescale of the first planetary 
encounter, t e , on the planet mass, M p , number of planets, 
N p , and the spacing between them, A m . When all planets 
have the same mass, they found that log 10 of the timescale 
is proportional to the planet spacing in Hill radii. When 
the number of planets N p > 5, the timescale is not strongly 
dependent on the number of planets. Here p = M p /M* is 
the planet mass in units of the stellar mass. The total mass 
in planets in a system would be M Pttota i — M P N P . 

For the number of planets N p = 3, the time to first en- 
counter, t e , varies from 10 4 to 10 6 years for planet masses 
fi ranging from 10 -5 — 10 -9 where a year represents the 
orbital period of the innermost planet. Large variations 
in planet mass are required to significantly change the 
time scale required for insta bility to develop (e.g., see Figure 
4 bv lChambers et al.ll 19961 ). The time to first planet/planet 
encounter is much more sensitive to planetary spacing than 
it is to planet mass since the hill radius is only dependent on 
mass ratio to the one third. This implies that a higher total 
mass in planets can be resident in a stable system when the 
mean planet mass is higher even though the planets must 
be somewhat further apart. A nearly unstable system with 
a larger number of Earth mass planets must have spacings 
such that the system has lower total mass in planets than 
a nearly unstable system containing small number of some- 
what more widely spaced Jupiter mass planets. For example, 
a planetary system that has not suffered a near encounter in 
10 6 years has spacing greater than distance in mutual Hill 
radii, A m > 7.5 for planet mass ratio fi — 10 -9 but only 
A m > 6 for fi = 10" 5 and . 

For a system observed with a dusty disk containing an 
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interior clearing, we assume that bodies residing between the 
planets must have been removed via dynamical processes 
on a timescale shorter than the age of the system. If nu- 
merous planetesimals lie interior to the clearing they could 
efficiently produce dust in the clearing. Dust produced in the 
outer disk can scatter into the clearing because the collision 
timescales estimated from the dust opacity is significantly 
shorter than the ages of the systems. Dust can also drift 
into the clearing because of radiative forces. Studies of dust 
particle integrations suggest that massive planets, at least 
Neptune in size, are req uired to account for a clearing in 
the dust distribution (e.g..lOzernov et al.l2000l:lLiou fe Zookl 
19991 ; iMoro-Martin fe Renull2002l ; iDeller fe Maddisonll200i 



Quillenll2006l . 12007ft . S harp eccentric edg es such as that in 
the Fomalhaut system ([Kalas et al.ll2005h i mply that clear - 
ing by planets is required in some systems (Quillen 2006). 
Our assumption is consistent with scenarios for the evolu- 
tion of our solar system. Most orbits in between Jupiter 
and Neptune becom e planet orbit cross ing within the age 
of the solar system (|Duncan et al.lll989T l. Nevertheless cra- 
tering history suggests that there was an epoch of early 
bombardment possibl y associated with the clearing of inter- 
planetary debris ( e.g. JComes et al.ll2005l ; lBottke et al.ll2005l ; 
IStrom et al.l2005l ) . We note that some systems with detected 
IR excesses may contain evidence f or inner exo-zodiacal belts 
(e.g., HD12039; [rimes et al.l 120061 ) in their spectral energy 
distribution or they may be present but below the detection 
limit. Nevertheless if planetesimals are originally widely dis- 
tributed and contain enough mass to form planets (as sug- 
gested by the cratering record in our solar system), then 
the lack of bright or detected exo-zodiacal belts implies that 
most planetesimals have been ejected, accreted onto plane- 
tary bodies or have fallen into the central star. 

In this paper we combine the assumptions that planetes- 
imals within a dust disk edge must be removed within the 
age of the system with estimates of the timescale t e for the 
first planet/planet encounter in a multiple planet system. 
While other explanations for clearing exist, such as a single 
planet migrating outward constantl y eroding away the inner 
edge of the disk (jGomes et al.ll2005f ). we only consider clear- 
ing via multiple planet induced instability. We estimate a 
minimum spacing between planets such that objects placed 
between them are likely to be removed on a timescale shorter 
than the age of the system while the planets themselves re- 
main stable. For simulations with planets, only the timescale 
until first planet/planet encounter depends on the planetary 
spacing and planet mass. Here we estimate the planet spac- 
ing for the instability of inter planetary objects by requiring 
the planets to be twice as far away as required for a first 
planet/planet encounter to occur during the age of the sys- 
tem. Such a system will have planets which will be stable 
on a timescale much longer than the age of that system. 
However, particles between the planet s will have timescale s 
shorter than the age of the system (e.g. JPuncan et aDll989ft . 
We expect that interplanetary particles residing in the sys- 
tem should have instability timescales similar to those es- 
timated from N-body simulations of planets spaced at half 
the spacing. By doubling the spacing, we ensure that un- 
der these timescales, the planets will be stable while the 
planetesimals will be cleared. The age of the system is a 
timescale for planetesimal clearing, not planetary instabil- 
ity. As all massive bodies influence each other and low mass 



bodies do not influence massive ones, we expect that our 
estimate for the required planetary spacing is conservative. 
Our spacing estimate will be larger than the actual one re- 
quired for efficient clearing of interplanetary bodies and so 
will lead to a lower limit on the number of, and total mass, 
in planets required to clear particles from clearings. 



2 NUMERICAL INTEGRATIONS 

Integrations were done using John Chambers' MERCURY 
package version 6. Mercury6 contains several N-body al- 
gorithms. We use the hybrid symplectic/Burlisch integrator 
since it is relatively fast and has the ability to compute close 
encounters. 

Ten massive bodie s of all the same mass ratio were used. 
IChambers et al.l (|l996T ) finds that as long as the number of 
planets in a system exceeds five, then the mean stability is 
nearly independent of number of planets. That is, a system 
of five planets has similar close encounter timescale depen- 
dence on separation as a system of ten or twenty planets. As 
we desire an estimate for the number of planets, we compute 
the stability timescale for a larger number of bodies in such 
a way as to be insensitive to actual number of planets. Each 
body has zero initial eccentricity, inclination, longitude, pe- 
riapse and mean anomaly. All planets were assumed to have 
density 1 g cm' 3 . They felt no non-gravitational forces. 

The initial semi-major axis, a n +i, of the n+ 1-th planet 
was initially set to depend only on the semi-major axis of 
the planet just interior or a n with a spacing 5 such that 

a.n+i = (l + 5)a„ (1) 

i + a(-) 

where th is the Hill radius for the n-th planet and A is 
the spacing in Hill radii. The spacings chosen for our inte- 
grations ranged from A = 2.5 to 11.0. This definition for 
pl anet spacing is somewh at different than the spacing used 
bv lChambers et all (|l996i ) . They use multiples of the mutual 
Hill radius rather than the Hill Radius. For small values of 
/x, the definitions are the same except for a factor of 2 1 ' 3 . 

We carried out a series of integrations for each planet 
mass ratio with different initial values of S setting the plan- 
etary spacings. Timescales for the simulation are given in 
orbital periods of the innermost planet. Integrations were 
calculated for ten million years or until a close encounter 
occurred. A close encounter is defined as when the distance 
between two planets was less than one Hill radius. When 
this circumstance occurred, the system was assumed unsta- 
ble and the integration terminated. The integrations did not 
include collisions. After every 1000 years, the six orbital el- 
ements for each planet was output and recorded. Planets 
were ejected from the system if their semi-major axis ex- 
ceeded 200 AU. Mass ratios of 10" 3 , 10~ 5 and 10~ 7 were 
used. We have used larger mass ratios than lChambers et al.l 
l|l996ft so we can consider systems with massive planets. 

The timescale to first close encounter time for our inte- 
grations are plotted as a function of planet spacing <5 / /i 1 ^ 4 in 

1 Available to download at: |http:/ /www. arm. ac.uk/~jecy] 
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log(t e ) = A + Blog(u/10" 7 ) + C(8/n' 



(1/4), 




Figure 1. The relationship between planet spacing, <5, (see equa- 
tion[2]l and time to first close planet/planet encounter, for systems 
with 10 coplanar planets. The spacings are multiplied by ^i 1 / 4 so 
that the sets of integrations with the same planet mass lie on lines 
with the same slope. Three sets of integrations are shown, those 
with planet mass ratio fi = 10~ 3 (pluses), 10 — 5 (x's) and 10 — 7 
(open squares). A fit to the three sets of integrations is shown 
with the solid lines. The fit is described with a single equation 
that relates timescale of first encounter to planet mass ratio and 
spacing (equation [2} . 



Figure HI IChambers et all (|l996t ) found a linear relationship 
between timescale of first encounter and planetary spacing 
in Hill radii but the slope of this relation is dependent on 
planet mass ratio. Only when rescaled by /x 1/4 (see their fig- 
ure 4) does the slope become independent of mass ratio. We 
confirm this here. 

We have fit the points shown in Figure[T]measured from 
our integrations with a single equation; 

log 10 (-^-) = -0.1 - log 10 W10- 7 ) + 3.2(5/ M 1/4 ). (2) 
yrs 

This relationship is shown with solid lines for the three mass 
ratios shown in Figure [1] The slop es of the lines in Figure [J 
are consistent with those found by [Chambers et all (ll996T ). 

The above equation is useful as it can be used to esti- 
mate a minimum spacing required for instability as a func- 
tion of system age. Rearranging equation [2] we find 



5 = 



3.2 



1 °Sio(rf:) + - 1 + lo gioWio~ 7 ) 

yrs 



(3) 



We confirm that S is not strongly dependent on planet mass, 
as we discussed in Section 1. Faults in this equation are 
discussed in the final section of this paper. 



3 ESTIMATE OF A LOWER LIMIT ON TOTAL 
MASS IN PLANETS 

We use the relationship measured from the integrations 
(equation[3]l to estimate the minimum spacing between plan- 
ets required for bodies to be removed within a clearing. We 
set t e in the above equation to the age of the system, re- 
quiring particle removal in the clearing within the lifetime 
of the star. Assuming that all planets have the same mass 



ratio, and that they are all spaced at twice 8 (defined by 
equation[2)|, the total number of planets, N p , must be larger 
than that satisfying 

Rout 



(1 + 2S) Np 



(4) 



where Rout is the semi-major axis of the outermost planet 
and Ri n is the semi-major axis of the innermost planet. 
The factor of two is used to ensure that over this timescale, 
the planetesimals are cleared while the planets themselves 
remain stable. Equation [3] for N p can be solved given 
Rin, Rout, and as 5 depends on t e and given the age of 
the star and an assumed planet mass ratio. 

The outermost planet must reside within the dust disk 
clearing so we can set R ou t to the radius corresponding to 
the dust temperature estimated from infrared spectra. We 
conservatively estimate Ri„ to be the distance at which wa- 
ter freezes (the iceline) based on the luminosity of the central 
star and given by: 



Rin — 



4-KcrTf 



(1/2) 



(•>) 



where T/ = 273K is the freezing point of water, L* is the 
luminosity of the star, and a is the Stefan-Boltzmann con- 
stant. The iceline represents a convenient limit for the inner 
radius. While massive planet formation models favor forma- 
tion outside the iceline, planets could migrate closer to the 
star. Massive planets could reside within the iceline so our 
estimate for the total number of planets is a lower limit. 
The number of planets in a system is then given by: 



[°gl()(Rout/ Rin) 

log 10 (l + 2J) 



(6) 



The factor of 2 for 8 on the right hand side follows because 
we have set the planets to be twice as far away as that 
required for instability in the pure 10-body system (see dis- 
cussion in section 1). 

The remaining quantity that we must choose to es- 
timate the number of planets is the planet mass ratio, 
/x. Recent studies suggest that maintenance of a clear- 
ing in the dust distri bution requires at l east a Nep- 
tune mass planet (e.g., Ozernov et al.l 2000l; Liou fc Zookl 



19991: iMoro-Martin fc Renull2002! : iDeller fc Maddisonll2005l : 
Quillenll2006l . 120071 ). We adopt this mass ratio, /i = 5x 10" 5 
as a starting point and then will consider the effect of varying 
this ratio. Once we estimate the number of planets present, 
the total mass in planets 



M p ,total > /J,M„N p 



( 7 ) 



The actual mass in planets is likely to be higher than the 
above estimate for a number of reasons: There may be mas- 
sive planets within the iceline. We have assumed that all 
planets have the same mass. If more massive planets are 
present then a larger total planet mass could be present. 
This follows because <5 is not strongly dependent on planet 
mass. Our integrations were run with ten planets but we 
may estimate that the number of required planets is fewer 
and systems with fewer planets are likely to be more stable. 

We apply our framework for estimating the minimum 
number and total mass in plan et s to the disks studied by 
IChen et al.l i|2006l ). IChen et all (|2006t ) has measured the 
radii of the clearings in the dust based on dust temperatures 
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measured from spectra observed with the Infrared Spectro- 
graph on board the SSC. Their sample consists of main se- 
quence stars which are located within 150pc and stars that 
are part of OB associations or moving groups. The proper- 
ties of these stars are listed in Table [1] with our calculated 
limits on the number of planets, N p , and total mass in plan- 
ets, M Pt totai within the clearings. We find that 3-6 planets 
of Neptune mass ratio are required and the total mass in 
planets is likely to be larger than ~ Q.2Mj where Mj is a 
Jupiter mass. 

Our estimates for the number of planets and total mass 
in planets depends on the assumed planet mass ratio, /j,. We 
investigate how the estimates depend on the assumed ratio. 
From the same stellar sample we have computed the mean 
number of planets using fi ranging from 10 -3 to 10 -7 . The 
results are plotted in Figure [5] The vertical line on this plot 
shows the location for Neptune's mass ratio. We find that 
the total number of planets is insensitive to the assumed 
mass ratio for fi > 1CP 5 . Planets with ratios lower than 
this are unlikely to be able to effectively eject material from 
central clearings within the lifetime of the system. While the 
total number of estimated planets is likely to be low, the 
total mass in planets could be an order of magnitude higher 
than given in Table 1. Jupiter mass planets with similar 
number and only slightly larger spacings would be capable 
ejecting material from the clearings during the lifetime of 
these systems. 

It is interesting to compute the total number of plan- 
ets for our solar system assuming a dust belt in the K uiper 
Belt were detected (as predicted bv lLiou fc Zo"okHl999h . Us- 
ing Rout = 30 AU and and an iceline at 1.0 AU for Ri n 
we would predict N p = 6.6 Neptune mass planets, giving a 
total planetary mass of 0.36 Mj. Much lower than the ac- 
tual planetary mass in our solar system in this semi-major 
axis range (1.4 Mj), we reiterate that our caclulations rep- 
resent an absolute lower limit on the total planetary mass 
of a system. Our solar system easily increases its planetary 
mass since the planets closer than Neptune have varying 
mass ratios larger than that of Neptune. Since stabilty is 
less dependent on mass ratio than spacing, fewer more mas- 
sive planets exist, thus increasing total planetary mass in 
the solar system. 



4 CONCLUSIONS AND DISCUSSION 

In this paper we have made the assumption that bodies, both 
planetesimal and dust particles, within a dust disk clearing 
must be removed from the clearing within the age of the 
system via gravitational interactions with planets. This as- 
sumption is consistent with some scenarios for the early evo- 
lution of our solar system. This assumption is combined with 
a rough estimate for the timescale for instability in multi- 
ple planet systems to estimate the planet spacing required 
for instability. The result is an estimate for the number of 
planets likely to reside in debris disk systems that host disks 
with clearings. We find that the number of required planets 
is 3-6, between the dust disk clear ing radius and th e iceline, 
for the sample of disks studied by (|Chen et al ] |2006h . assum- 
ing a planet mass ratio like that of Neptune. We find that 
the number of estimated planets is only weakly dependent 
on the assumed planet mass. At least 0.2 Mj total mass 
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Figure 2. We consider the average number of planets for the 
systems in Table [T] for assumed mass ratios between 10 — 3 and 
10 — 7 . The curved line, oc (log 10 /i) -0 25 , is a best fit line for 
the average number of planets in a system based on mass ratio. 
We find that the number of estimated planets is insensitive to 
the assumed planet mass ratio for planets sufficiently massive to 
effectively empty a clearing. If more massive planets are present 
then the total mass in planets could be an order of magnitude 
higher than listed in Table [T] 

in planets is required in each system. The planets could be 
more massive and an order of magnitude more mass in plan- 
ets could reside in these systems. 

Our estimate neglects planets that lie within the ice- 
line but planets could migrate or form closer to the star. 
The estimate for the planet spacing is based on integra- 
tion of a 10-body system but our estimated planet number 
falls below 10. Further numerical experiments show that 
if the number of planets is decreased to five, a decrease 
in S/fi 1 ' 14 ' by ten percent results in the same instability 
timescale. The estimate for spacing is twice that for insta- 
bility in a 10-body system but really we require an instabil- 
ity timescale for m assless bodies in betw een massive bodies 
(e.g., as studied bv lDuncan et al1ll989h . We have assumed 
that the planet/planet encounter timescale is related to a 
clearing timescale and this is not necessarily true as small, 
low mass objects do not feel forces on each other the way 
that large, massive bodies do. A more sophisticated treat- 
ment would require integration of a much larger number of 
particles and low mass as well as planetary mass bodies. 
We could also examine the effects of resonances on stabil- 
ity timescales and the phase space of systems with massive 
planets as was done bv lChatteriee et al.l l|2007l ). Such an ex- 
ploration could improve on the crude constraints on multi- 
planetary systems based on observations of debris disks that 
we have explored here. Here we have assumed that central 
clearing of debris and dust is primarily due to scattering 
by planets. However collisional cascade models for dust pro- 
duction in debris disks predict that the dust production rate 
drops with time and that the cen ters of disks evolve faste r 
than than their outer parts (e.g., iDominik fc Decml 120031 ) . 
It is possible that future observations and statistical studies 
may differentiate between collisional cascade evolution and 
planetary clearing models. 
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Table 1. Stars with Dusty Disks with Clearings 
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The stellar masses (M*), outer disk radii (R ou t), Luminosity (L*), and stellar ages are 
given for the sample disc ussed and studied by IChen et al ] (l2006ll. "Age of HR 5314 3 from 
IZuckerman fc Sond ||2004T| . ''A ge of HR 95086, HP 146897 from ldeZeeuw et alj dl99Sh . Note 
as can be seen from Table 1 of lChen et al ] l|2006h there are discrepancies between available 
age estimates. The rightmost two columns are calculated using Equations \E\ and [3] and show 
the estimated number of planets residing between the iceline and the clearing edge (at Rout)- 
The iceline radius Ri n is estimated using equation \E\ Each planet is assumed to have the 
mass ratio of Neptune. The number of planets, JV P , is estimated assuming that the planets 
are sufficiently close together that all interplanetary debris has suffered a close encounter 
with a planet within the lifetime of the system. The total mass in planets M p tota i = N p fiM* 
is a lower limit for the total mass in planets residing in the system. Our own solar system is 
included to emphasize this point. As is the case in our Solar System, if more massive planets 
are present, then more mass could be present in planets between the dust disk and the iceline 
even though the planets must be further apart. 



